Hydrologic Design 4501
Chapter 6: Surface Water and Runoff Processes
Ardeshir Ebtehaj
University of Minnesota
1- Watersheds and Runoff
A waterhsed, a catchment or a drainage basin, are three synonymous terms that define a locus of points on the Earth's surface that drain precipitation water to a single point called the wathershed outlet.
Left: Schematic of precipitation in a watershed flowing the precipitation water towards the outlet point. Right: A schematic of how watersheds are nested and watershed size varies depending on the basin outlet selected (right, credit: Marsh, 1998, p. 170)
2- Basics of Surface Runoff
During a rainfall event, there are two types of water storages:
- Detention Storages: Short-term storages that are depleted by overland flows and infiltration.
- Retention Storages: Long-term storages that are depleted by evaporation.
Surface runoff or overland flow is defined as the portion of rainfall, snowmelt, and/or irrigation water that runs over the soil surface toward the stream network rather than infiltrating into the soil. The definition of runoff also may include the interflow and/or return flow, which together with the surface runoff make up the volume of stormwater or total runoff. The return flow is referred to the shallow subsurface lateral infiltrated flow of water that returns back to the surface.
During a precipitation event, the detention storages begin to fill out and water starts to flow over land surfaces. Overland flows join together over a hillslope and form more concentrated flows that channelize the underlying soil creating the so-called channel or stream network. Several channel networks eventually joint and drain into the main stream and generate the streamflow.
A schematic of the main runoff components in the absence (left) or presence (right) of snowmelt.
The flow that you see in a stream or river typically consists of three main hydrologic fluxes that we briefly covered in Chapter 1. These fluxes are: overland flow, return flow, and baseflow as shown in Fig. 2. The overland flow component can be split into two categories:
- Hortonian overland flow: The overland flow that occurs when precipitation rate is greater than infiltration capacity. It is also called infiltration excess runoff.
- Saturated overland flow: The overland flow that occurs when soil becomes saturated from below due to rise in the groundwater table. The produced runoff through this mechanism is also called saturation excess runoff.
Infiltration exceeds precipitation (left), Hortonian overland flow (middle) and Saturated overland flow (right).
As we discussed, the infiltration rate is the water flux that enters the soil at the surface. It is often expressed as depth of water per time, for example, 10 millimeters per hour.
If precipitation rate is much higher than the infiltration capacity, overland flow occurs immediately after the onset of precipitation. Infiltration excess is commonly observed with short-duration intense rainfall. It also occurs most often over land surfaces with high clay content or where the surface has been altered by soil compaction, urbanization, or fire.
However, when precipitation rate is less than or equal to infiltration capacity, no overland surface runoff occurs. For prolonged precipitation event, the water table may rise and produce the saturated overland flow from below. It is most common with long-duration, light-to-moderate rainfall, or because of successive precipitation and or snowmelt events at upstream. Saturation excess runoff can occur anywhere, when the soil is saturated from below. It is the most common runoff generation mechanism in humid climates with gently sloped or flat basins. The saturate overland flow and return flow often occur near the stream channels, where the groundwater depth is often shallower and it rises quickly in response to infiltration -- called groundwater ridges.
A schematic of overland versus saturation excess runoff over the hillslope and near the channel.
3- Streamflow Hydrograph
The streamflow hydrograph is the time series of flow rate at a specific location on a stream. These measurements are typically captured at stream gauges, which measure flow rate based on river stage and a flow rating curve. The USGS is the primary source of stream gauge data in U.S.
For simplicity, we will categorize two different hydrographs typically analyzed based on the time scale of interest. An annual hydrograph focuses on the flow over a time period of a year at daily time steps (shown below). A storm hydrograph focuses on the flow corresponding to a precipitation event for time periods ranging from hours to weeks depending on the watershed size and storm duration.
An example of an annual (left) and storm (right) hydrograph for Vermillion River near Empire, MN. (Credit: USGS)
There are multiple characteristics of a basin and its river network that determine the shape of the outflow hydrograph including:
- (a) size of the drainage basin
- (b) slope
- (c) hydraulic roughness
- (d) natural and channel storages
- (e) stream length
- (f) channel density
- (g) antecedent soil moisture
- (h) other factors such as land cover
The effects of basin characteristics on the flood hydrograph (Masch 1984)
The dynamics of precipitation events also affects the shape of the hydrograph.
Effects of the shape of the rainfall hyetograph, storm size and movement on the shape of the streamflow hydrograph.
3-1 Hydrograph Components
Hydrographs are typically divided into two components for analysis. The baseflow, which is the long-term contribution of groundwater flow to a stream hydrograph and the direct runoff, which is the portion of the hydrograph that is in response to an excess rainfall of an event and is also referred to as quick flow.
The storm hydrograph can be divided into a few segments as follows:
Main components of a storm hydrograph (McKerchar 1993). Note that the saturated overland flow and interflow (return flow) are the main components of the direct runoff hydrograph.
3-2 Baseflow Separation
Three different methods for basflow separation are:
- Straight line method: Draw a line from the beginning of the rising limb to end of the falling limb, where the recession starts.
Components of the streamflow hydrograph: (1-2) baseflow recession, (2-3) rising limb , (3-5) crest segment, (4) peak flow, (5-6) falling limb, and (6-7) baseflow recession. The straight line method is used for baseflow septation. At the beginning of the falling limb is an inflection point (5), where the direction of the curvature changes.
- Fixed-base method: Draw a tangent line from the pre-event recession curve to the peak time and then connect it to a point on the hydrograph that the surface runoff is assumed to be negligible.
Fixed-based method for baseflow septation.
- Variable-slope method: Draw a tangent line from the pre-event recession curve to the peak time and a tangent line from the post-event recession curve to the inflection point of the falling limb. A straightline is drawn between the peak discharge and inflection point.
Variable slope method for baseflow septation.
3-3 Excess Rainfall and Direct Runoff
The direct runoff determined from a hydrograph does not exactly correspond to the volume of precipitation delivered over the basin since there are losses to storages. Therefore, when predicting direct runoff from precipitation inputs, we divide the rainfall hayetograph into:
- Excess Rainfall: The amount of rainfall that is neither retained on the land surface nor percolate into the groundwater as a recharge flow. This rainfall volume flows across the watershed as surface runoff or return flow to generate direct runoff. The excess rainfall is also known as the effective rainfall.
- Abstraction: The rainfall that is absorbed primarily by infiltration followed by groundwater recharge as well as interception by vegetation and surface storages. The excess rainfall hyetograph (ERH) is equal to the total rainfall minus total abstraction.
The conceptual relationships of rainfall, infiltration rate, and excess rainfall are shown in the following figures.
The concept of rainfall excess, which is the difference between the actual rainfall hyetograph and losses largely due to the infiltration. Generally, we either assume a constant rate (red line) or a time varying infiltration (black) rate in computation of the excess rainfall hayetograph (ERH).
Approximating excess or effective rainfall is crucial for stomrwater managment and prediction of streamflow for future rainfall events. There are different methods for determining effective rainfall depending on whether streamflow data is available or not.
3-3-1 Excess Rainfall using Streamflow Data
To determine the excess rainfall hyetograph (ERH), the first step is to separate the baseflow from the direct runoff hydrograph (DRH). One of the most common methods for obtaining the ERH is the ϕ-Index method. This method assumes a constant rate of abstraction throughout the storm event and is calculated as follows:
where
is the excess runoff depth,
is the observed rainfall depth over a time interval
with a decreasing order
, where M is the number of rainfall pulses that contribute to direct runoff, and
is called the constant abstraction rate.
A schematic representation of the ϕ-index method. The
and thus
are the results of direct runoff, which is obtained by subtracting the baseflow from the total observed flow.
The process consists of checking the above equation for descending values of
until a positive ϕ value is achieved. Then double check that the depth of excess rainfall equals the depth of direct runoff. --------------------------------------------------------------------------------------------------------------
Example Problem 6.1: We have rge following information of a rainfall storm over a watershed area of 7.03
with a baseflow of 400 [cfs]:
Solution:
is obtained by subtracting the baseflow 400
from the streemflow (column 3) data resulting the data in column 6Total amount of excess rainfall uniformly distributed over the watershed is
Let us begin the iterative solution:
not feasible
not feasible
feasible!
Note: If
then we had to try
and continue the process.
For characterizing the amount of excess rainfall from stream data, we can also calculate the runoff coefficient (C):
From the previous example, we can easily calculate the runoff coefficient as:
--------------------------------------------------------------------------------------------------------------
The runoff coefficient is used extensively in engineering practices to determine the peak flow rate of a storm event based on average rainfall intensity over small watersheds through the rational method:
where
is the peak flow rate [cfs], i is the design rainfall intensity
, and A is the basin area
. Note that the dimensions for the rational method do not match and it is purely an empirical relationship. In practice, this method is used to design the size of catch basins and storm sewers by calculating the peak flow rate for a certain return period of the design rainfall event. There are standard estimates of the runoff coefficients based on different soil types, land use and slopes. For example for lawns (C=0.05-0.35), forest (0.05-0.25), concrete streets (0.7-0.95).
--------------------------------------------------------------------------------------------------------------
Example Problem 6.2: Consider a design problem where the peak discharge is required to size a storm drain inlet for a 2.4-acre parking area where the design rainfall intensity is 8.6 in/hr and the runoff coefficient of is 0.95.
Solution:
--------------------------------------------------------------------------------------------------------------
Where drainage area is characterized by distinct subareas that can be characterized by different runoff coefficients, a weighted runoff coefficients should be determined. The weights are based on the area of each land use and is formed by the following equation.
in which
is the area for land cover j,
is the runoff coefficient for area j, nis the number of distinct land covers within the watershed and
is the weighted mean runoff coefficient. Therofere, in this case the peak runoff is
.
--------------------------------------------------------------------------------------------------------------
Example Problem 6.3: Compute the runoff for the following land use chracteristics in heterogenous basin for design rainfall of 3.6 
Solution:
--------------------------------------------------------------------------------------------------------------
3-3-2 Excess Rainfall Using Infiltration Methods
If streamflow data is not known, then we must resort to infiltration models to calculate the excess rainfall hyetograph. It is common to use infiltration models like Green-Ampt to estimate the ERH. The key assumption is that Direct Runoff = Precipitation - Infiltration.
Recall that the Green-Ampt infiltration rate is defined as:
where
denotes changes of infiltration rate in time t and
is the cumulative infiltration from 0 to time t. Additionally, we showed that the cumulative infiltration can be computed as follows after the ponding time:
Now for the purpose of creating the ERH, we must define cumulative infiltration after the ponding time. To that end, we showed that:
Where the cumultive infiltration at the time of ponding (
), was defined as follows:
As described in the following, finding ERH using infiltration method is an iterative process.
Algorithm for determination of infiltration and ponding time under variable rainfall intensity:
WHILE
(
: the rainfall duration)
If infilteratoion rate is greater than precipitation (
) (no ponding at time t)
Else ponding occurres, when infiltration rate
is smaller than precipitation rate
(
) .
End
END
--------------------------------------------------------------------------------------------------------------
Example Problem 6.4: A rainfall hyetograph is given in col(1) and (2) of the following table. The soil is sandy loam with initial saturation
,
and
. Determine the excess rainfall hyetograph (ERH).
Solution:
@ t=0 [min]
@ t=10 [min]
The same process and calculatoin continue with no ponding, as shown in the following excel sheet, until
@ t=60 [min]
We need to use the Green-Ampt formula for computation of
:
Cumulative excess rainfall = cumulative rainfall - cumulative infiltration = 2.41(col 3) -2.21(col 6) = 0.2 [cm](col 7)
After 140 [min], ponding ceases as
An excel file of the above calculation table is available in the class website.
Below is a plot showing the results of the example problem. Notice that all abstractions before the time of ponding (where surface runoff is generated) is labeled initial abstractions. After
, everything above the green line is effective rainfall and everything below is infiltrated (losses). Solution plot to Green-Ampt example.
3-3-3 SCS Method for Runoff Calculation:
The old USDA Soil Conservation Service (SCS; now called USDA-NRCS) developed a runoff estimation method based on the following hypothesis.
Schematic of the hydrograph components for the SCS method.
P: Total precipitation [m, in]
: Initial abstraction [m, in]
: Excess precipitation [m, in]
: Water retained in the watershed [m, in]S: Maximum retention storage [m, in]
The above equation can be written as follows:
and thus
From data of field experiments for small watersheds, we can assume
, which leads to the following expression of excess precipitation:
Experimental results allowed parameterization of S using a curve number, which describes the retention capacity of differing landscapes:
[inches]where CN is called the curve number 
CN=100 (Impervious surfaces, e.g., pavements)
CN<100 (Natural surfaces)
Solution of the SCS runoff equations for normal soil moisture condition (from U.S. Department of Agricultural Soil Conservation Service 1972)
The CN depends on:
- The antecedent moisture condition (AMC)
- Soil type
- Land use
The curve numbers in the above figure or available tables are for normal antecedent moisture conditions (AMC II). The soil moisture conditions for SCS method are defined as follows:
Definitions of AMC based on 5-day antecedent rainfall
Therefore, correction factors are necessary for the dry and wet conditions as follows:
The CN values also depend on the soil type with four categories:
- Group A: Deep sand, deep loess, aggregated silt (lower CN values, higher retention)
- Group B: Shallow loess, sandy loam
- Group C: Clay, loam, clayey loam
- Group D: Soils with heavy content of clay (high CN values, low retention)
Runoff curve numbers (average watershed condition,
, USDA, 1986)
The curve number in the above tables are for average size watersheds (
) and normal antecedent moisture condition (AMC II). For watersheds with several sub-catchments with different CNs, the area-averaged composite values of CN shall be computed. --------------------------------------------------------------------------------------------------------------
Example Problem 6.5: Compute runoff from 5 inches of rainfall on a 1000-acre watershed. Hydrologic soil type is 50% (B) and 50%(C). The watershed land use is:
- 80% residential area (40% soil type B, 40% soil type C) that is 30% impervious
- 20% paved roads (10% soil type B, 10% soil type C) with curbs and storm sewers
Solution:
for wet condition:
and then we an compute S and
for this new curve number to account for anrecedent soil wetness.
5- Runoff Hydraulics
Surface runoff in a watershed occurs first as a thin sheet of overland flow in the upper slopes for a short distance (< 100 ft) and then eventually produce or merge to the channel flow that you see in gullies, streams and rivers. The goal for this section is to briefly quantify the properties of overland sheet flow.
Overland flow to the streams and modeling concepts
Here p is precipitation intensity
or
, θ is the slope angle,
is infiltration rate
or
, V is average velocity
, L is slope length
, and y is the flow depth
. Let's write a mass balance or continuity equation for the above control volume as follows:
- Inflow:
(precipitation) - Outflow 1:
(infiltration) - Outflow 2:
(discharge per unit width)
Conservation of mass (inflow = sum of outflows) results in
=
.
Note that in the above expression, the unit is in terms of discharge per unit width (e.g.,
). Thus, the overland flow per unit width at the end of the slope is:
Recall that from the Darcy-Weisbach equation, we have
- f: friction coefficient
- L : pipe length
- D : pipe diameter
- V: velocity
: head loss.
By definition, the hydraulic radius in a pipe or a channel is defined as
A: wetted area [
P: wetted perimeter [
For example, for a rectangular channel we have
(when
). Recall that for a pipe
and thus
, when the the Darcy-Weisbach equation is used for a rectangular channel
5-1 Laminar Overland Flow
Next we need to use the conservation of momentum to define the average velocity and height of the overland sheet flow. These equations depend on whether flow is laminar (parallel streamlines) or turbulent (mixed by eddies), which can be determined by the flow Reynold's number:
(pipe)
(wide rectangular channel)
where ν is kinematic viscosity
. We know from fluid mechanics that generally the flow is laminar when
. Recall that from the Moody diagram, we have
for laminar flow in pipes; however, for uniform laminar overland sheet flow, the following formula shall be used (Chow 1998):
where p is the precipitation rate [
].
Recall that the Energy Grade Line (EGL) is
and for a uniform flow, we have
, where
is the slope of the EGL,
is the head loss over the horizontal length L.
From the Darcy-Weisbach equation, for a uniform flow,
and
then we have
where
.
5-2 Turbulent Overland Flow
Once overland flow becomes turbulent (
2000), the roughness factorf becomes independent of the Reynolds number, and the Manning's equation shall be used to explain the velocity of overland sheet flow as follows:
where n is Manning's roughness coefficient and
is the hydraulic radius for a wide open channel. The Manning coefficient, is an empirically derived coefficient, which is dependent on many factors, including surface roughness and channel sinuosity. Substitute
into Manning's equation one can obtain:
Typical values for the Manning's coefficients n for overland flow are given on the next slide. These values are valid for sheet flow with depths less than 0.1 ft. We will revisit the Manning equation in the next chapters.
Typical values of the Manning roughness coefficient (n) for overland sheet flow.
Hydraulic radius formulas for simple open channel geometries. These formulas will be revisited in the next chapters.
5-3 Travel Time
One main reason for defining the velocity of these runoff mechanisms is to determine the travel time in a basin. As we know, for different segments of the river with the same velocity, the travel time can be approximated as follows:
Every point within a watershed will have a travel time based on its flowpath to the outlet. When rainfall begins, the areas with the shortest travel times are contributing first to the outflow. As the storm progresses, larger areas of the watershed contribute to the outflow until the time when all areas in the watershed are contributing to flow at the outlet. This time scale is known as the time of concentration (
) and can also be defined as the travel time of precipitation water from the longest flow path to the outlet.
Schematic showing the growing contributing areas and the time of concentration.
Approximate average velocities in
of runoff flow for calculating time of concentration. The unconcentrated condition addresses the flow velocity in the upper part of the watershed prior to the overland flows accumulating in a channel, while concentrated condition refers to the flow velocity in a channel network (From drainage manual, Texas highway department, 1970).
Example: Calculate the time of concentration of a watershed in which, the longest flow path covers 100
of pasture at slope 5% and a 1000
long rectangular channel having 2
width with Manning coefficient
and slope of the chennel is
%. The channel receives a lateral flow of 0.0096
. The flow velocity is 3
for the 100
of pasture with slope 5% (see the above Table) and thus
.
For the rectangular channel with length 1000
, we have
The total discharge for segment
is
and
. Thus, after finding y from the Manning equation, we can compute
. We solve the problem for river segments
in the following table: For example, for the first segement
. The depth of the flow in third row of the table is obtained by solving for the root of the Manning equation for a rectangular channel
. The depth are calculated for the flow rate in the second row for
ft. After calculating the depth, we caculate the velocity as
. The
represents the average velocity over each stream interval. For example,
for the second segment.